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Abstract 

Based on the Lue-Starkman conjecture on the dynamical screening of the brane 
cosmological constant in the DGP scenario, we extend this proposal to a general DGP- 
inspired F(R, <f>) Model. We show that modification of the induced gravity and its 
coupling to a quintessence field localized on the brane, affects the screening of the 
brane cosmological constant and also phantom-like behavior on the brane. We extend 
our study to possible modification of the induced gravity on the brane and for clari- 
fication some specific examples are presented. As a result, phantom-like behavior can 
be realized in this setup without violating the null energy condition at least in some 
subspaces of the model parameter space. The key result of our study is the fact that 
a DGP-inspired F(R, <j)) scenario has the best fit with LCDM and recent observations 
than other alternative theories. 
PACS: 04.50.-h, 98.80.-k 
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1 Introduction 



Recent evidences from supernova searches data [1,2], cosmic microwave background (CMB) 
results [3-5] and also Wilkinson Microwave Anisotropy Probe (WMAP) data [6,7], show an 
positively accelerating phase of the cosmic expansion today and this feature shows that the 
simple picture of the universe consisting of the pressureless fluid is not enough to describe 
the cosmological dynamics. In this regard, the universe may contain some sort of the addi- 
tional negative-pressure dark energy. Analysis of the WMAP data [8-10] shows that there is 
no indication for any significant deviations from Gaussianity and adiabaticity of the CMB 
power spectrum and therefore suggests that the universe is spatially flat to within the lim- 
its of observational accuracy. Further, the combined analysis of the WMAP data with the 
supernova Legacy survey (SNLS) [8], constrains the equation of state Wd e , corresponding to 
almost 74% contribution of dark energy in the currently accelerating universe, to be very 
close to that of the cosmological constant value. In this respect, a LCDM ( Cosmologi- 
cal constant plus Cold Dark Matter) model has maximum agreement with the recent data. 
Moreover, observations appear to favor a dark energy equation of state, Wd e < —1 [H]- 
Therefore, a viable cosmological model should admit a dynamical equation of state that 
might have crossed the value Wd e = —1 in the recent epoch of cosmological evolution [12]. 
In fact, to explain positively accelerated expansion of the universe, there are two alternative 
approaches: incorporating an additional cosmological component ( dark energy) in matter 
sector of the general theory of relativity ( G^ = 8nG(T^ + T^ ark ^) where and 
Tj£ ark ^ are energy-momentum tensor of ordinary matter and dark energy respectively), or 
modifying geometric sector of the theory (dark geometry) (G^ + G^ rk = 8nGT^) at the 
cosmological scales. Multi-component dark energy with at least one non-canonical phan- 
tom field is a possible candidate of the first alternative. This viewpoint has been studied 
extensively in the literature ( see [13,14] and references therein ). Another alternative to 
explain current accelerated expansion of the universe is extension of the general relativity to 
more general theories on cosmological scales. In this view point, modified Einstein-Hilbert 
action via f(R) -gravity ( see [15] and references therein) or braneworld gravity [16-18] are 
studied extensively. In this framework the geometric part of the Einstein's field equations 
are modified. For instance, DGP ( Dvali-Gabadadze-Porrati) braneworld scenario as an IR 
modification of the general relativity explains accelerated expansion of the universe in its 
self-accelerating branch via leakage of gravity to extra dimension. In this model, equation 
of state parameter of dark energy never crosses the uj(z) = —1 line, and universe eventually 
turns out to be de Sitter phase. Nevertheless, in this setup if we use a single scalar field 
(ordinary or phantom) on the brane, we can show that equation of state parameter of dark 
energy can cross phantom divide line [19]. One important consequence in the quintessence 
model of dark energy is the fact that a single minimally coupled scalar field has not the ca- 
pability to explain crossing of the phantom divide line, = — 1 [20]. However, a single but 
non-minimally coupled scalar field is enough to cross the phantom divide line by its equation 
of state parameter [13,14]. Lorentz invariance violating vector fields in an interactive basis 
are other possibility to realize cosmological line crossing [21]. Lue and Starkman [22] based 
on the analysis firstly reported by Sahni and Shtanov [23] have shown that one can realize 
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the phantom-like effect ( increasing of the effective dark energy density with cosmic time) 
in the normal branch of the DGP cosmological solution without introducing any phantom 
field. This type of the analysis then has been extended by several authors [24] . The normal 
branch of the model which cannot explain the self-acceleration, has the key property that 
brane is extrinsically curved so that shortcuts through the bulk allow gravity to screen the 
effects of the brane energy-momentum contents at Hubble parameters of the order of the 
inverse of crossover distance [22]. Since in this case H(t) is a decreasing function of the 
cosmic time, the effective dark energy component is increasing with time and therefore we 
observe a phantom-like behavior without introducing any phantom matter. It is important 
to note that crossing of the phantom divide line in this viewpoint is impossible without 
introduction of a quintessence field on the brane [24] . This idea has been studied further to 
incorporate curvature effects [25]. The importance of this type of reasoning lies in the fact 
that we don't need to introduce phantom fields that violate the null energy condition and 
suffer from several theoretical problems. 

Here we are going to study phantom-like effect in the normal branch of a general DGP 
inspired F(R, 0) scenario. The DGP inspired F(R, <fi) scenarios have been studied in Refs. 
[26,27]. Our motivation to study phantom-like behavior of this extension of the DGP scenario 
is the fact that to have crossing of the phantom divide line on the DGP brane we have to 
incorporate a quintessence field on the brane [24]. On the other hand, it is reasonable to 
assume that induced gravity on the brane can be modified. In fact, as has been argued in 
Refs [28], generalized version of DGP scenario ( such as modified induced gravity), can be 
ghost free and can give rise to transient acceleration ( see also [23] and [29]). Here we are 
focus on the normal branch of the scenario which is ghost-free. We show that for the case 
with F(R,<j>) = |(1 - i<p 2 )R, the effective dark energy density reduces by increasing the 
values of the non- minimal coupling, £. We extend our study to the general f(R) -gravity 
to explore the role played by the modification of the induced gravity on the screening of 
the brane cosmological constant and the phantom-like effect. We show that phantom-like 
behavior can be realized in this setup without violating the null energy condition at least in 
some subspaces of the model parameter space. The key result of our study is the fact that 
a DGP-inspired F(R, 0) scenario has the best fit with LCDM and recent observations. 



2 Non- minimal DGP Cosmology 



2.1 The Setup 

The action of the DGP scenario in the presence of a non-minimally coupled scalar field on 
the brane can be written as follows [27] 

/rrr^ f / 777,^ I \ 

d 5 x^^TZ+ J d*xV=q{-fa{4>)R[q\- -<TV '^W ^ - V (cf)) + m 3 4 K + C m j , 

L - - J y=0 

(1) 

where we have included a general non- minimal coupling a((f>) in the brane part of the 
action( for an interesting discussion on the possible schemes to incorporate NMC in the 
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formulation of the scalar-tensor gravity see [30,26], and also [31] for a braneworld viewpoint). 
m\ = (87rG) _1 and y is the coordinate of the fifth dimension and we assume that brane is 
located at y — 0. qab is five dimensional bulk metric with Ricci scalar TZ, while is 
induced metric on the brane with induced Ricci scalar R. K is trace of the mean extrinsic 
curvature of the brane defined as 



Kf„, = ^ lim ( 



o V 



K 



+ 



K 



\1V 



y=+e 



(2) 



and corresponding term in the action is York- Gibbons- Hawking term [33] (see also [34]). 
The ordinary matter part of the action is shown by the Lagrangian C m = C^q^, ip) + 
where ip is matter field and corresponding energy-momentum tensor is 



T, 



-2 



.5Crn 

5q^ 



mi 



(3) 



and A is the brane cosmological constant. Note that we assume that in addition to brane 
cosmological constant, there is some quintessence scalar field localized on the brane to have 
a more general framework and in order to realize phantom divide line crossing. The pure 
scalar field Lagrangian, = — \q^ v V v 4>— V(<j>), yields the following energy-momentum 
tensor 

V = V M 0V,0 - -<V(V0) 2 - q^V{<t>). (4) 
The Bulk-brane Einstein's equations calculated from action (1) are given by 

m\ (k A b ~ -^QabT^ + 



a{4>) (iV - \% V R) ~ V M V,a(0) + g M ,D (4) a(0) 5{y) = 8 A li 8 B v T llv 8{y), (5) 



where is 4-dimensional (brane) d'Alembertian and T^ u 
be rewritten as follows 



T/ju, + t uv . This relation can 



ml (Kab ~ \g AB n) + mla(<j>)5 A *6B v (r^ ~ \q^R) <%) = 8 A ^ B v T^{y) (6) 
where T^ v is the total energy-momentum on the brane defined as follows 

% v = mlVpVMfi ~ mlq^n^a^) + T^, (7) 

From (6) we find 



Gab = TIab — -g A B^ 



and 



G^u — I R 



\%vR) 



X, 



mla((j>) 



(8) 
(9) 
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for bulk and brane respectively. The corresponding junction conditions relating the extrinsic 
curvature to the energy-momentum tensor of the brane, have the following form 



lim 



K 



y=+e 
y=-c 



1 



mla(4>) 



ml 



Rnv — -q[Mvq a/3 Ra/3 



(10) 



Now we study cosmological dynamics in this setup. Since DGP scenario accounts for embed- 
ding of the FRW cosmology at any distance scale [33,34], we start with following line-element 



ds 2 = q^ v dx^dx u + b 2 (y, t)dy 2 = —n 2 (y, t)dt 2 + a 2 (y, t)^ijdx l dx^ + b 2 (y, t)dy 2 . 
In this relation 7^ is a maximally symmetric 3-dimensional metric defined as 



^ij &ij ~t~ k 



X 2 X j 

1 — kr 2 



'ID 



12) 



where k = —1,0, 1 parameterizes the spatial curvature and r 2 = XiX 1 . By computing com- 
ponents of Einstein's tensor and using junction condition given in equation (10), we arrive 
at the following Friedmann equation in this non-minimal DGP setup [27] 



H 2 + X = 



a 2 3m 2 a((f)) 



p m + p<p + P\ + Po 



2 S 

Pm + P<f, + PA -m§Q!(0)-j 

a 



1+eJl + — 
V Po 



• (13) 



where po 



6m| 



Pa 



= m\a(4>)h and p rn is density of ordinary matter on the 



brane. Also, e — ±1 shows the possibility of existence of two different branches of DGP- 
inspired FRW equation corresponding to two different embedding of the brane in the bulk. 
Neglecting the dark radiation term |f ( where £ is an integration constant) which decays 
very fast at late-times, we rewrite equation (13) as follows 



H 2 = 



8nG 



, \ A 1 /ll rSvrG/ x A 

( Pm + ^) + 3 + 2^ + £ te + ^^( Pm + ^) + 3 



(14) 



where r is a crossover distance defined as r = ^_dgp Q; (0) an d G = G e // — l/87rm|o;(0). 
Now, the Friedmann equation (14) can be rewritten as follows 



H 2 = 



8ttG 



A 



f \ 1Y H 



3 \ r rv > 3 
We use this equation in our forthcoming arguments. 



(15) 



2.2 Lue-Starkman Screening of the Brane Cosmological Constant 

As we have pointed out in the introduction, Lue and Starkman have shown that one can re- 
alize phantom-like effect, that is, increasing of the effective dark energy density with cosmic 
time, in the normal branch of the DGP cosmological solution without introducing any phan- 
tom field. The normal branch of the model which cannot explain the self-acceleration, has 
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the key property that brane is extrinsically curved so that shortcuts through the bulk allow 
gravity to screen the effects of the brane energy-momentum contents at Hubble parameters 
H ~ Tq 1 where r$ is the crossover distance [22]. Since in this case H{t) is a decreasing 
function of the cosmic time, the effective dark energy component is increasing with time and 
therefore we observe a phantom-like behavior without introducing any phantom matter that 
violate null energy condition and suffers from several theoretical problems. In the first step, 
in this section we study the phantom-like effect in the normal branch of a DGP inspired 
non-minimal scenario. In other words, here we suppose that there is a quintessence field 
non-minimally coupled to the induced gravity on the DGP brane. We emphasize that we 
have included a canonical (quintessence) scalar field to incorporate possible coupling of the 
gravity and scalar degrees of freedom on the brane. This provides a wider parameter space 
with capability to handle the problem more complete. In fact, inclusion of this field brings the 
theory to realize crossing of the phantom divide line [24]. As has been shown by Chimento et 
al. , the normal branch of the DGP scenario has the capability to describe phantom-like effect 
but it cannot realize crossing of the phantom divide line without introducing a quintessence 
scalar field on the brane. With this motivation, here we have considered the existence a 
canonical scalar field on the brane that couples non-minimally with induced gravity. In the 
next section we incorporate possible modification of the induced gravity on the brane too. 
Considering the normal branch of the equation (15) with e = —1, we have 

rr2 SttG/ x A H 

H 2 ^—( Pm + P A + ---. (16) 

6 v ' 6 Tq 

Comparing this equation with the following Friedmann equation 

^= 8 -f{p m + P,) + 8 -f P m\ (17) 

87rG Jeff) _ A H / 1R x 

— pDE - J- To (18) 

Existence of a quintessence field nonminimally coupled to the induced gravity on the brane 
leads to a redefinition of the crossover scale as r = £dgp Q; (0)- Using definition of r , 
equation (18) can be rewritten as follows 

^ G Jeff) _ A a-'H 

-^~Pde - t - - B • U y J 

<J <J tDGP 

Now we assume a conformal coupling of the scalar field and induced gravity as follows 

a{ct>) = \{l-if). (20) 



we find 1 



1 Note that this comparison is not perfect since G in equation (16) is an effective quantity defined as 
G = G e ff = 1/8^77?! a (</>). However, since screening of the brane cosmological constant can be attributed 
just to the last two terms of the right hand side of equation (16), this comparison is actually possible. 
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The values of the £ is constraint by the observations from different viewpoints ( see for 
instance [35,36]). The division by 1 — £0 2 in our field equations unavoidably introduces the 
two critical values of the scalar field ±0 C = ±— 7=, for £ > 0, which are barriers that the 

scalar field cannot cross. Note that in these values, the effective gravitational coupling, its 
gradient, and the total stress-energy tensor diverge ( see [31] for more details). 
Now, by adopting ansatz (20), equation (19) can be rewritten as follows 



8ttG 



(eff) 



-Pde 



A 
3" 



2H 



'DGP{ 



£0 2 



(21) 



Figure 1 shows the variation of versus £ in a constant time slice. In plotting this 

figure we have used the ansatz a{t) = a t u , <fi(t) = 0o^ _M with v = 1.2 ( an accelerating 
phase of expansion) and /i = 0.9 ( a decreasing quintessence field). The range of £ are 
chosen from [36] constraint by the recent observations. As this figure shows, by increasing 
the values of the nonminimal coupling, p%^]P decreases in a fixed time slice. Figure 2 shows 



Pot 




Figure 1: Variation of p D J versus £ in a t = constant slice. 

the variation of the effective dark energy density versus the cosmic time. As this figure 
shows, pjfg increases with cosmic time and this is exactly the phantom-like behavior we 
are interested in. Note that this phantom-like effects is realized without introducing any 
phantom matter on the brane and only screening of the brane cosmological constant causes 
such an intriguing effect. Although the existence of a canonical scalar field non-minimally 
coupled to the induced gravity on the brane has no considerable effect on the phantom-like 
behavior but as figures 1 and 3 show, increasing the values of the non-minimal coupling leads 
to the reduction of the effective dark energy on a constant time slice. Figure 3 shows the 
variation of the effective dark energy versus the cosmic time and non-minimal coupling. We 
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Pde 5 ~ 




-10 J 

Figure 2: Variation of the effective dark energy density versus the cosmic time. The effective 
dark energy density increases with time and therefore shows a phantom-like behavior. 

note that while the introduction of a phantom field requires the violation of the null energy 
condition, here this energy condition is respected since we have not included any phantom 
matter on the brane. Since the phantom-like dynamics realized in this setup is gravitational 
( the quintessence field introduced here plays the role of standard matter on the brane), the 
null energy condition cannot be violated in this case. 




t o.o 



Figure 3: Variation of the effective dark energy versus the cosmic time and the non- minimal 
coupling. For a constant £ the model realizes phantom-like effect versus cosmic time. 
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3 DGP-inspired F(R, (f)) Gravity 



3.1 The Setup 

Now we extend our previous analysis to the more general case with DGP-inspired F(R, (f>) 
models. In other words, we incorporate possible modification of the induced gravity on the 
brane. We assume also a general coupling between a quintessence field localized on the 
brane and modified induced gravity ( these types of theories have been studied extensively 
and from various perspectives, see for instance [26, 30, 37]). The action of this model is as 
follows 

3 2 i 

S=^ Jd 5 (x)^^ + Jd\x)^(^F(R, ( f>)--q^V, ( f>VA-V( ( f>)+mlK + L m )m 

where the first term shows the usual Einstein-Hilbert action in 5D bulk with 5D metric de- 
noted by qab and Ricci scalar denoted by 9ft. The second term on the right is a generalization 
of the Einstein-Hilbert action induced on the brane. This is an extension of the scalar-tensor 
theories in one side and a generalization of /(i?) -gravity on the other side. We call this model 
as DGP-inspired F(R, (f>) scenario, y is the coordinate of the fifth dimension and we suppose 
that brane is located at y — . q^ is induced metric on the brane which is connected to qab 

A 
SttG 



via q^ v = 5^ A 5u B gAB- We denote matter field Lagrangian by L m = L m (q fJ , 1/: 'ip) + with 



energy-momentum tensor defined as = — 2|^j + q fll/ L m . The pure scalar field lagrangian 
is L<f, = — \q^ v V ^(pV v 4> — V{4>) which gives the following energy- momentum tensor 

V = V M 0V,0 - -^(V0) 2 - q^V{<t>). (23) 
The field equations resulting from this action are given as follows 

"' ' ~ \gAB^) +m\5 A V(R^ V - \q^R)5(y) = 5 A ^b "(f^+T^+T^)6(y). 



F'(R,< 

( 24 ) 

In this relation T^ u = p,^^ where T^ v is the energy-momentum tensor in matter frame and 

f MI/ = pTjjffi ■ A prime denotes differentiation with respect to R. Also, TjfJ 17 ""^ is defined as 
follows 



T^ rv) = pJ^^ [l^ ( F (R, 0) - RF'(R, 0)) + (V(i?, 0)) ( q m q v p - q^q af3 
In the bulk, T A b = and therefore 



(25) 



and on the brane we have 



Gab = ^ab - \g AB ^ = (26) 



1 q; 

Gnu = Rav — T.QlJ.vR = o ' (27) 

2 f«3 
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where T^ v = f^ v + + T^ arv \ The corresponding junction conditions relating quantities 
on the brane are as follows 



lim \K 



y= 



ml 



1 



^F'{RA) 
m% 



al3 R, 



aj3 



y=0 



(28) 



A detailed study of weak field limit of this scenario within a harmonic gauge on the lon- 
gitudinal coordinates and using Green's method to find gravitational potential, leads us to 
a modified (effective) cross-over distance in this set-up as follows ( see [27] for details of a 
similar argument) 



m 



IfdF 



2m 3 d \dR 



(29) 



where as usual £dgp = 



3.2 Cosmological Implications of the Model 

As we have explained in the previous section, embedding of FRW cosmology in DGP setup 
is possible in the sense that this model accounts for cosmological equations of motion at 
any distance scale on the brane with any function of the Ricci scalar. To study cosmology 
of a DGP-inspired F(R, (f>) scenario, we consider the line element as defined in equation 
(11). Also, we assume that the scalar field <fi depends only on the cosmic time on the 
brane. Choosing a Gaussian normal coordinate system so that b 2 (y,t) = 1, non-vanishing 
components of the Einstein's tensor in the bulk plus junction conditions on the brane defined 
as 

" (dF\ / a cr na 
^dR'^n 2 a 2n 2 a 2 n 3 a 2a 2 



hm Q [d y n} y y : 



2nm\ 



ml 



y=0 



n 



3ml 



>dF\ 
^~dR< 



lim [^ a] jK(0 = 4 [(^) (-¥- + ") 
:^+o L y iy mll^dR J ^n 2 a nJ 



(tot)^ 










y= 


i 










■/dF\p^a 




- y=0 




XdRJ 3ml _ 


y=o 



(30) 



(31) 



yield the following generalization of the Friedmann equation for cosmological dynamics on 
the brane ( see [26,27] for machinery of calculations for a simple case) 



H 2 + 



k 



a 2 3m 2 3 F'{R, 0) 



(f> tot + Po 




2 




1 + 6a 


1 + — 


ptot 






Po 





m 2 3 F'(R, <f))eo 



(32) 
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where e = ±1 shows two different embedding of the brane, po = 



6m|_ 
miF'{R4 



and Eq = 



3^ — a' 2 + kja 2 is a constant with respect to y ( with a' = ^) ( see [34,27,38] for more 
detailed discussion on the constancy of this quantity). Total energy density and pressure are 
defined as p^ to£) = p + p^ + p curv + pa and p^ tot "> = p + p^ + p CM ™ + p A respectively. The 
ordinary matter on the brane has a perfect fluid form with energy density p and pressure p, 
while the energy density and pressure corresponding to non-minimally coupled quintessence 
scalar field and also those related to curvature are given as follows 



p v = 



\^ + nV(0) - 6^//0 



(33) 



also 



V<t> 



1 • 



2 dF /V n : , dFH ■ 2 d 2 F ., 

■(f) 



2 - \/(0) + - — (0 - -0) + 4— —<j> + 



2n 2 n 2 d0 n dd> n 2 n 2 dd) 2 



y=o 



(34) 



P 



(curv) 



mi, (\ 



F'(i?,0)V2 



F(R,(/))-RF'(R,<f)) -3RHF"(R, (/>)), 



(35) 



P 



(c«to) 



/? ■ [2RHF"(R, 0) + HF"(i2, 0) + i? 2 F'"(i?, 0) - - [F(i2, 0) - i*F'(i2, 0) ) .(36) 



where = tt£4t is the Hubble parameter on the brane. Ricci scalar on the brane is given 



by 



a(0,t) 



R — 3— H — - 



a 2 n 2 



a v a 



a\2 an' 
j -6 — 

a n. 



Note that cosmological dynamics on the brane is given by setting n(0,t) = 1. With this 
gauge condition we recover the usual time on the brane via transformation t = /*n(0, rj) drj 
where 77 is conformal time. It is interesting to note that the equation of state parameter of 
the scalar field defined as 



2RH + R ] F"(R, if>) + R 2 F"'(R, 



±F(R, 4,) + iSF'(R, 4>) 



\<f 2 + V(4>) - 6 



crosses the phantom-divide line u; = — 1 in the favor of recent observations [26,30]. 



(37) 



3.3 The Phantom- Like Behavior 

Now in this DGP-inspired F(R, 0) model, the crossover scale takes the following form 

_<F%R 1 0) _ 
f ~ 2m* ~ ^ ' ' 



(38) 
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also 



Po 



3m 2 F'(R, 0) 



2i) 



Neglecting the dark radiation term in equation (32), we find 



H 2 



SnG(p m + P4> ) A J_ 

3 + 3 + 21) ^ " 



4 + «2 



1 r 87rG( Pr , 



+ 



A 



(39) 



where G = G e // = (8-KmlF' (R, 0)) and a prime denotes differentiation with respect to R. 
By adopting the negative sign we find 



__ 2 87rG(p m + P<p ) A ^%R>) 

-n — ^ r — ; 



H 



DGP 



(40) 



We can compare this equation with equation (17) to conclude that the screening effect on 
the cosmological constant is modified as follows 



8nG ff _ A 
~z~Pde — ~W 



F'(R, 



H 



I 



DGP 



(41) 



As an important especial case, for F(R, 0) = f(R) we find the screening effect in a general 
/(-R)-gravity 



8ttG 



Pde — A 



f'(R) 



H 



3 zdgp 
Now as an enlightening example, we set for instance 

1 



(42) 



F(R, 



'i - Z4> )[R - (i - n)C{R/C) n ] , 



where ( is a suitably chosen parameter ( see for instance [15] and [26]). With this choice, 
one recovers the general relativity if n — 1. For n ^ 1, we obtain from equation (40) 



8ttG 



eff 

Pde 



A 



2H 



3 rua 3 e DGP (l-C(l) 2 )(l-n(l-n)C^ 1 -^R n - 1 ) 
For spatially flat FRW geometry the Riici scalar is given by 

R = 6- + 6(- f. 

a a 



(43) 



(44) 



To have an intuition of phantom-like behavior in this case, we adopt a suitable ansatz so that 
a(t) = a t u and <fi(t) = 0o^ _/1 - We set v = 1.2 and /i = 0.9 that are reliable from physical 
grounds. Figure 4 shows the variation of p®^j versus n in this DGP-inspired F(R, 0) model. 
As we see, phantom-like behavior can be realized for n > 0.73 and n < —0.60. In other 
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Figure 4: Variation of the effective dark energy versus n in a DGP-inspired F(R, (f>) model with 
F{R,4>) = |(1 - ^ 2 )[R - (1 - n)( 2 (R/( 2 ) n }. Phantom-like behavior can be realized for n > 0.73 
and n < —0.60. 
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Figure 5: Variation of the effective dark energy versus n and the non-minimal coupling 
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words, for —0.6 < n < 0.73 the effective dark energy in this DGP-inspired F(R, 0) model 
has no phantom-like behavior. 

Figure 5 shows the variation of the effective dark energy versus n and the non-minimal 
coupling. By increasing the values of £, the effective dark energy density reduces but for a 
fixed value of £, there is phantom-like effect for appropriate values of n. Also, figure 6 shows 
the variation of the effective dark energy versus n and the cosmic time. The phantom-like 
effect ( increasing the values of the effective dark energy) can be realized for suitable range 
of n. 



rho(DE) 



-11- 




Figure 6: Variation of the effective dark energy versus n and the cosmic time. 



3.4 The Expansion History 



To investigate expansion history of our model and comparing it with other alternative theo- 
ries, we study luminosity distance versus redshift in this scenario. For a dark energy model 
with constant equation of state parameter, the luminosity distance versus redshift can be 
expressed as follows 



:i+z) 



H^dz 



+ zf + (1 - + z)^+") 



nl/2 



(45) 



and for a LCDM model, u = — 1. Now, the evolution of the cosmic expansion in our DGP- 
inspired F(R, 0) model is given by 

H(z) 



1 



1 



+ 4^o 



M 



+ 4fi0[(l + 2)3(1+0,) _ 1] 



, (46) 
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where by definition tq = £dgpR'{R,4>)- The luminosity distance versus redshift in a LDGP 
model can be expressed as [22] 

'VM-WL'WY (47) 

and in our DGP-inspired F(R, <fi) scenario, this quantity denoted as d'[ DGP (z) is given by 

^- w = (1+2) /;^, (48, 

where H(z) is given by equation (46). Figure 7 shows a comparison between expansion 
histories of LCDM, LDGP and our FDGP scenario for F(R, <j>) = |(1 - £(f) 2 )[R - (1 - 
n)( 2 (R/( 2 ) n ] with £ = 1/6 and n = 0.8. Note that this value of n lies in the appropriate 
range required for realization of the phantom-like effect obtained in the previous subsection 
and it is also suitable for describing late-time acceleration ( see for instance the paper by 
Sotiriou and Faraoni in Ref. [15]). A LCDM scenario has very good agreement with recent 
observations. As we see here, the FDGP scenario is closer to LCDM more than LDGP. 
In other words, FDGP has better agreement with recent observation than LDGP. Therefore 
FDGP provides a better framework for treating phantom-like cosmology without introducing 
any phantom field. Since we have not introduced any phantom matter on the brane ( (f) is 
a quintessence field which plays the role of standard matter on the brane), it seems that 
the null energy condition should be respected in this setup. However, as we will show in 
subsection 3.6, this is valid only for some specific values of the model parameters and only 
in some subspaces of the model parameter space. 

LDGP 
- . - . - LCDM 
FDGP 



4 n 




Figure 7: Luminosity distance versus redshift for three alternative scenarios. For a model universe 
with F(R,(j)) = |(1 - £<f> 2 )[R - (1 - n)C 2 (R/C 2 ) n ] with f = 1/6 and n = 0.8, there is better fit 
between FDGP and LCDM. 
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3.5 Dynamics of the Equation of State Parameter 



To have more detailed discussion on the cosmological dynamics in this model, we find from 
equation (40) the following relation ( with 87rG = 1) 



H 



--^-iM*^*))- 



(49) 



Considering the energy conservation equation which is expressed here as p tot + 3H(p tot + 
Ptot) = where p tot = p m + P<p and p tot = p m + p , we find 



H = — 
2 



p tot (l +Utat) 



e DGP [F>(R,< 



l + 4[F'(R,<p)}% GP (^f + l) 



|. (50) 



There is no superacceleration in this DGP-inspired F(R, <fi) scenario if the following condition 
holds 

P <ot (l + - tot )> w[m0)]2 - (51) 

To have a general relativistic interpretation of the expansion history of this model, we rewrite 
the energy conservation equation as follows 



p ef f + 3H(l+UJ e ff)p e ff = 

and using equation (40) we have 



(52) 



Peff 



F'(R, 0) 



3H 

+ 



-2 . 



F'(R,(f>) F'(R,(f>) 



^DGP d-DGP 

By comparison of equations (52) and (53), we find 



(53) 



1 +u, 



eff 



F'(R, 0) 



-l • 
H 



-2 . 



F'(R,<j>) F'{R,<j>) 



H^DGPPeff 



^DGPPeff 



(54) 



To realize the phantom phase in this DGP-inspired F(R : 0) model, the condition l+u e ff < 
should be fulfilled. This leads us to the following condition: 



H < F'(R,, 



H F'{R,<j>) 



(55) 



It is obvious that this model has the potential to describe the crossing of the phantom divide 
line. 
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3.6 The Null Energy Condition 



It is important to check the validity of the null energy condition in this setup. In fact, the 
main feature of this setup is the realization of the phantom-like behavior without introducing 
any phantom matter on the brane. The null energy condition is respected if the condition 
p + p > is valid. In our case, this condition is given by p tot +Ptot > where p tot and p tot are 
defined in the subsection 3.2. Figure 8 shows the variation of y = {p + p)tot versus n for some 
specific values of redshift. As this figure shows, there are appropriate subspaces of the model 
parameter space that the null energy condition is respected in this setup. This is enough 
to say that this DGP-inspired F(R, <fi) model realizes the phantom-like behavior without 
violating the null energy condition, at least in some subspaces of the model parameter space. 
For instance, at z = 0.25 ( which is corresponding to the epoch of the phantom-divide line 
crossing), the null energy condition is respected if n < 1.8. Albeit, those values of n are 
adequate that are supported observationally( by, for instance, solar system tests). It should 
however be noticed that this range seems more restrained at higher redshifts. The reason 
for violation of the null energy condition in some subspaces of the model parameter space 
lies in the fact that a modified theory of gravity of the form f(R) is equivalent to a theory of 
standard gravity plus a scalar field. With f(R) gravity, we have shown that one can mimic 
a phantom-like behavior without introduction of a phantom field, but when the scenario is 
written in the Einstein frame, the resulting scalar will violate the null energy condition. So, 
it is natural to accept that in our model there are some subspaces of the model parameter 
space that the null energy condition can be violated. The main achievement is however the 
existence of other subspaces that respect the null energy condition. 

z=l 

z=0.25 

z=0 



lOO-i 




y 50- 






^i- 


-2 -1 0- 
-50- 

-100- 


i • ^\ 3 

\ \ \ 
\ \ 
\ \ 

I 



Figure 8: y = (p + p)tot versus n for some specific values of redshift. There are some subspaces 
of the model parameter space that null energy condition is fulfilled for a model universe with 
F(R, </>) = i(l - ^ 2 )[R - (1 - n)C 2 (R/C 2 ) n ] with £ = 1/6. 
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4 Summary 



Based on the Lue-Starkman conjecture on the dynamical screening of the brane cosmolog- 
ical constant in DGP scenario, in this paper we have extended this proposal to a general 
DGP-inspired F(R, 0) Model. Firstly, we have studied phantom-like behavior in the normal 
branch of an extension of DGP model where a quintessence field is coupled non-minimally to 
the induced gravity on the brane. The reason for incorporation of this canonical scalar field 
lies in the fact that without scalar field it is impossible to realize phantom divide line crossing 
in DGP setup. We have shown that the effective dark energy density decreases by increas- 
ing the values of the conformal coupling £ in a constant cosmic time slice. However, for a 
constant £, we have phantom-like behavior ( increasing of the effective dark energy density 
with cosmic time) in the normal branch of the scenario without introducing any phantom 
field. Then we have extended our study to a general DGP-inspired F(R, 0) scenario where 
we incorporate possible modification of the induced gravity on the brane. In this case we 
obtained some new and interesting results which we summarize as follows: by adopting the 
ansatz F(R, <fi) = \{l — i4> 2 )[R — (1 — n)( 2 (R/( 2 ) n ] , we have shown that phantom-like be- 
havior can be realized in the normal branch of the scenario if n > 0.73 and n < —0.60. In 
other words, for —0.6 < n < 0.73 the effective dark energy in this DGP-inspired F(R, 0) 
model has no phantom-like behavior. Investigation of the expansion history of this model 
shows that this DGP-inspired F(R, 4>) scenario has the best fit with the recent observational 
data. In fact this model is very close to a LCDM scenario. Finally we found conditions 
for transition to phantom phase of this model which has the potential to realize phantom 
divide line crossing. For the case of a quintessence scalar field non-minimally coupled to the 
induced gravity on the brane, the null energy condition is fulfilled since there is no phantom 
matter on the brane and the phantom dynamics is essentially gravitational which saves the 
null energy condition. Also that the brane tension does not violate the null energy condition 
too. For a general DGP-inspired F(R, <fi) scenario, the null energy condition is respected 
only in some subspaces of the model parameter space depending on the choice of the model 
of modified gravity. 
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